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The electromagnetic field around a Kerr black hole inside a current loop is sometimes used as the 
basis of a toy model for discussing the properties of particle orbits near astrophysical black holes. 
The motivation for the present paper is to correct the pubhshed solution to Maxwell's equations 
with a charged current loop. Dipole approximations and closed-form expressions in the extreme 
Kerr hmit are also presented. Using the corrected solution, it turns out that imposing a vanishing 
electromotive force produces a loop with a potential which is finite everywhere outside the black 
hole. Ring solutions can be combined into solutions with multiple rings or current discs. 



;3 

1—5 

oo 



^ I. INTRODUCTION 

Attempts have been made over the years to explain the generation of untrarelativistic particles through physical 
effects around black holes located inside accretion disks (see [I] for a review). In some situations, large electric 
C potentials can be produced by dynamo action if the black hole is rotating in a magnetic field generated in the 
accretion disk. In the original Blandford-Znajek model, for example, the induced electric fields drive currents and act 
as the source of energy fluxes [2^ i3j . 

In the context of some of these models, it becomes important to study particle orbits and particle interactions in 
electric and magnetic fields close to the black hole. These studies are increasingly being replaced by magnetohydro- 
O dynamics (MHD), but the particle picture is needed when a fluid description is not appropriate. The study of circular 
orbits, for example, gives the location for the inner edge of the accretion disk for free charges [IHB]. In the presence 
$-H of the magnetic field, the equations governing the charged particle orbits are no longer seperable, and the particle 
motion can be come quite complicated outside of the equatorial plane [Tj [8] . This may well have an effect on issues 
such as charge separation and the validity of fluid approximations. 
^Sj The large potential around the black hole implies that Penrose processes can be an important source of untrarela- 

^ tivistic particles 0. Familiar limits on the efficiency of the Penrose process [TUlfTT] are not valid for charged particles, 
and the particles emerging from the ergoregion can have energies close to the available electric potential p^. For 
^1;^ example, magnetic fields of 10^ Gauss around 10^ solar mass black holes can produce positrons and electrons with 
^—^ energies as high as lO^^eV through reactions such as 77 — ?► e+e~ inside the ergoregion. 

All of the discussions of particle motion mentioned above have used the field from a current loop, or its dipole limit. 



as a toy model to represent the magnetic field. The solution to Maxwell's equations for the uncharged current loop 
was obtained a long time ago by Chitre and Vishveshwara [T3] , and generalised to a charged current loop by Petterson 
^ |14j . The general class of stationary solutions was analysed by Bicak and Dvorak [T51, who also considered some 
• • special cases including current loops, dipoles and point charges. Surprisingly, the result of Petterson is not consistent 
^ ^ with the other two papers, which is unfortunate because this is the only published solution which solves for the vector 
k> potentials that are needed to study the particle dynamics. 

j_J In order to explain the differences further, consider the general vacuum axisymmetric solution to Maxwell's equa- 

te tions, with a set of four arbitrary constants for each angular mode. These constants can be related to the current and 
charge by matching vacuum solutions inside and outside the current loop. The inconsistency between the three papers 
mentioned above is due to an error in Ref |14| when matching the solutions. (There is also a minor typographical error 
in the constants of Ref. [H]). Correcting the values of these constants leads to physical consequences, for example the 
potential at the horizon of the black hole changes and the zero EMF condition used in some astrophysical scenarios 
occurs for a different value of the black hole charge. The errors in Ref. [M] have subsequently been reproduced in 
standard texts "T. The motivation for the present paper is to present the correct electromagnetic potentials for a 
charged current loop. 

The current loop solutions are appropriate when there are no free charges surrounding the black hole and there 
is no radial current flow. In real situations the region surrounding the black hole may contain particle-antiparticle 
pairs which conduct current freely. This type of model could be built up by superimposing radial currents with the 
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ring solutions. The simplest possibility is to have radial current flowing out along the equatorial plane and then back 
along the symmetry axis, as in the toy model by Li 16 . 

There are problems with the simple Blandford-Znajck model, for example it is likely that most magnetic field lines 
pass through the accretion disk TT^ and the disk is then subject to a variety of MHD instabilities [18 . On the other 
hand, a similar model may still hold in systems with relatively low accretion rates where the inner parts of the disk 
are effectively Keplerian orbits 19 . These issues are not considered further here, and the results below can be used 
when the electromagnetic fields generated by the accretion disk can be replaced by the fields generated by current 
loops. 



II. MAXWELL FIELDS IN THE KERR BACKGROUND 



The general solution to Maxwell's equations on a Kerr background is described in the textbook by Chandrasekhar 
[2D] . In principle, we ought to solve the full Einstein-Maxell equations, but the gravitational back- reaction of the 
magnetic field can be ignored when BM << 1 in gravitational units, or 

5 « 2 X 10^" G (1) 



in Gauss. In typical astrophysical situations, the electromagnetic field strengths are likely to be insufficiently large to 
have any significant effect and the Kerr metric for a rotating black hole can be used consistently. 
The Kerr metric in Boyer-Lindquist coordinates (t, r, 9, (p) is given by 

,9 A o sin^ 9 S , o „ ,„o 

ds^ = --ujI + —^ojI + -dr^ + Y.de^ (2) 

2-1 2-i Za 

where 

(jJi = dt — asinddcf), u)2 — [t"^ + a^)d(f) — adt (3) 
E = cos^ 9, A ^ r'^ + a'^ ~ 2mr (4) 

The Maxwell field strength i^^^ is represented by three complex scalar quantities 

00 = F^J^m'' (5) 
= Ff.^il^n'' + m^'*m-')/2 (6) 
cl>Q = F^,m''*n'' (7) 

where the Newman-Penrose tetrad vectors are 

= i(r2 + a2,A,0,a) (8) 
= ^(r2-Ha^-A,0,a) (9) 
= 3^ (iasin6', 0, 1, icosec6') (10) 

and p = r + i cosO. 

The stationary axisymmetric solutions to Maxwell's equations with vanishing sources are determined by solving the 
Teukolsky equation 

- A*,„ - cosec6'(sin 6'*,e),e + (cot^ + 1)-^ = (11) 

where commas denote derivatives and the field 4* = —A(f>o. The general solution for 0o can be expressed in terms of 
angular mode functions Pi^ {cos 9), 

^ CO 

(po = ^Y.i"'^P'Au)Pi\cos9) + PiQiAu)Pi\cos9)) (12) 

The coefficients a; and f3i are complex constants and P/(u) and Qi{u) are Legendre functions of u = {r — M)/{rj^—M)^ 
where r+ is the radial coordinate of the event horizon. The solution can be checked using the identity 

{APur{u))^r = l{l + l)Pi{u). (13) 
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The terms involving are regular at the horizon and the terms involving Qi{u) are regular at spatial infinity. 

Typically, we have different inner and outer solutions separated by a region with charges and currents. 

The reconstruction of the gauge potentials given in ref [20] is incomplete, but with a little extra work it is possible 
to extend the method used there to rederive the inner solution given by Petterson [Ti] 



At = - ^ 1 -^PiA'^)Pii^'^^^){^oii - a cos 0ai) + ^ Pi (u) sin 6 Pi^ {cos 9) {a cos 9 + r a}) \ 

1=1 ■' 

A,^fsinH + ±{ 
1=1 ^ 



(14) 



°^Pi^r{u)Pi{cos9)sm^{9){ral - acos9al) (15) 

+ ^ ^° P; (u) sin 6>P;^ (cos g) (a cos 6>a[ + - Pi,riu)Pi\cos9) sm9a{\ (16) 

/ (t + 1) J 

where Q is the black hole charge, at is a constant and ai — a] + ia[. The possibility of a constant term in A,f, is 
excluded by the regularity of the vector potential along the axis of rotation. 

The outer solution uses the alternate set of Legendre functions, but it is also important to account for the charge 
of the source region q. Since this will be important later, it is worth noting the full outer solution, 

At - - ''^^^ '^^ ~Y.{ ^QiA^)Pi (cos 9){r^l - a cos 9pl) + [u] sin gP;^ (cos 9) {a cos 9p\ + r/3?)| (17) 
1=1 ^ ^ 

- "^(^+ sin^ e + ^ |^Q;^,(^)P;(cosg) sin\9){rpl - acos9p]) (18) 
1=1 ^ 

+ ^^^^ Qi{u)siii9Pi\cos9){acos 9(31 + rfi\ ) - j^^-^ Qi.r{u)Pi\cos9)sin9l3\^ (19) 
The gauge has been chosen so that Ai vanishes at spatial infinity. 

III. CHARGED CURRENT LOOP 

The field due to the charged current loop is obtained by matching the inner and outer solutions at the position 
of the loop. The safest way to do this is to use the Maxwell scalars 0o or (/)2, because these are readily expressed 
in terms of orthogonal polynomials in cos 6*. Matching as was attempted in Ref il4J, leads to errors because 
contains both Pi{cos9) and Pi^(cos6'), and these are not orthogonal functions. 

The current density for the charged current loop in the Boyer-Lindquist coordinates is given by |15] 

J''= ( ;^(5(r-ro)(5(cos6i),0,0, 4^(''-^o)'5(cos6i) ) . (20) 







The parameter q is equal to the total charge of the loop, as defined by an integration of the charge density over a 
surface S of constant time, 

' J^dS^ = q. (21) 



Gauss' law, applied to a large sphere, implies that this is the same charge which appears in the potential (18). 

The relationship between the Maxwell scalars and the source is obtained through the corresponding Teukolsky 
equation, 

- A*,„ - cosec6'(sin6'*,9)^e + (cot^ 6* + 1)* = AttJ (22) 
The source function is related to the current density by 

J = S(A + /i)S~V,„. - Y.{5* + vr - t*)E-V„, (23) 
where A and 5* are directional derivatives along n and m* respectively, and 

A ia sin 9 ia sin 9 

u = , TT = =, r = (24) 
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For axisymmetric sources, we may express the fields in angular modes with components Ri, where 



Ri^ J 'i' Pi^ {cos e)d cos 9, 



and similarly for J;. These components satisfy 



A A 



In the source-free regions, (f>o is given by Eq. (12), and the integration (25) yields 

-A2l{l + 1) 



Ri = 



V2 21 + 1 
-A2l{l + 1) 
^21 + 1 



aiPi^r{u) inner region 
PiQi.r{u) outer region 



When the current density ( 20 ) is substituted into eq. ( 23 ) , the differential equation ( 26 ) for Ri becomes 

l{l + l) 



R 



i.rr 



A 



-Ri = F{r)6ir - ro),, + G(r),5(r - ro). 



where 



Fir) = _^4^/(!i±i^!)^p,i(0) 
V2 ^-0 



Integrating Eq. ( 29 ) through the delta function gives 



Combining these equation gives the solution for 



1 + 1/21 ^^^2GM , , V2F{rol ,1 



where w is the Wronskian, 



w = PiMQi^rM - Qi{uo)Pi^r{ua) 



r+- M 



The corresponding solution for a; is obtained by replacing Pi by Qi. The real and imaginary parts of /3; 
identified by using Eqs. (30) and (31), 



PI 



1 + 1/2 1 4TTaI-2q 
l{l + l)w ro 
1 + 1/2 1 ( 47r/(r2 + a^) - 2aq 
lil + l)w \ 



Similarly, for ai, 



a, 



roA{ro 
1 + 1/2 1 ATraI-2q 



Pi,rMPi{0) 

Piiuo) 



AttI 
l{l + l) 



PiAuo)\ Pi\Q) 



l{l + l)w ro 
1 + 1/2 1 UT:I{rl + a^) - 2aq 
l{l + l)w\ roA(ro) 



Ql,rMPl{0) 

Qiiuo) 



AttI 
1(1 + 1) 



QiAuo)}Pi\o) 



(25) 
(26) 

(27) 
(28) 

(29) 

(30) 
(31) 

(32) 
(33) 

(34) 

(35) 
can be 

(36) 
(37) 

(38) 
(39) 
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These results are consistent with Chitre and Vishveshwara in the q = hmit, allowing for the different definition 
of /. The real parts differ from Petterson [14j[2T] by a factor of tq, but the imaginary parts are totally different. 
The constant term at remains to be evaluated, and for this we use the field along the axis of rotation, 



At 



At 



A 



7-2 + a' 
A 

+ a' 



1=1 

oo 



r{Q + q) 



+ OLt 



(40) 
(41) 



1=1 



For at, we require continuity of At aX r — tq. From Eq. (32), we get 



Evaluating the sum and using Eq. ( 30 ) gives 



y^aF^ ^ / + 1/2 



2TTla 



at 



(42) 



(43) 



Note that the term at represents the value of the potential at the horizon of an uncharged black hole. 

The magnetic field lines (in any synchronous frame) can be traced by the contours of constant A^, and a couple 
of examples are shown in Fig. [T] The plot on the left shows the field lines in a vertical plane through the axis for 
an extreme Kerr black hole with a current loop at r = 6Af, and the plot on the right shows the field lines for an 
extreme Kerr-Newman black hole (with a small charge). The former case is an example where the black hole expels 
the magnetic flux at the horizon. This black hole version of the Meisner effect has been noted before [22l - f24] . The 
plots shown here demonstrate clearly that extreme black holes can have magnetic flux across the horizon when they 
are charged. In fact, A^ = Q/2 at the black hole equator and = at the north pole. 



z/M 




-10 



z/M 




-10 



8 10 



FIG. 1: The magnetic field lines are shown in a vertical plane through the axis for an extreme Kerr black hole with an equatorial 
current loop at r = 6M (left), and for an extreme Kerr-Newman black hole (right). The mode sum has been truncated at 
I = 16, and the small irregularities ar r = 6M are due to this truncation. 
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IV. CHARGE ACCRETION 

The physical interpretation of the horizon potential is that the rotation of the black hole combined with the magnetic 
field of the loop generates a large electromotive force (EMF). In an astrophysical context, as first pointed out by Wald 
[55] . free charges can migrate along the magnetic field lines to cancel the EMF, resulting in a charged black hole. 

If the total system of black hole plus current loop is to remain neutral, then charges should also migrate to the 
current loop. We shall refer to the resulting solution as the zero-EMF solution with vanishing total charge. The 
zero-EMF and vanishing total charge conditions where analysed in Ref. [T3|, but these results used an incorrect 
form for the electromagnetic potential. The charge on the current loop was corrected by Linet [57] , but without 
obtaining the vector potential. (Note that Linet uses a different definition of the current /.) The full solution is given 
below for the first time. 



A. Zero-EMF solution with vanishing total charge 



It will be sufficient to consider the electrostatic potential (|40j) along the axis of rotation. The gauge has been chosen 
so that the potential vanishes at infinity, and so the total EMF is given by the potential —At at the horizon r_|_. It 
vanishes when 



at 



Q_ 
2M' 



(44) 



The quantity at was evaluated in (43). If we require vanishing total charge then q = —Q and the loop has a charge 

inMIa 



q 



ro - 2M 

Examples of the vector potentials for the zero-EMF solution are shown in Fig. [2] 



(45) 



-0.2 



At 




-0.8- 



FIG. 2: The vector potential solutions in the equatorial plane for a loops with radius 6M around an extreme Kerr black hole. 
The dipole solution for a field of strength B is also shown. The fields are scaled by the modulus of the horizon value. 
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It is interesting to note that there is a lower hmit to the radius of the loop in the zero-EMF solution when a > M/2. 
Suppose that the loop contains particles of charge e and iV_ particles of charge — e, so that 

q = {N+-N_)e. (46) 

The current is bounded by having all of the charges moving at close to the speed of light, hence [50] 

{N+ + N^)e 



\I\ < 



27rro 



The charge q = —Q is related to the current by (45), and the limit (47) implies that 

ro > M +{2Ma- M'^fl'^ 



(47) 



(48) 



when a > M/2. This is larger than the radius of the minimal circular stable orbit for black holes close to the extremal 
limit. 



B. Finitenes of the potential 

It is natural for the electromagnetic fields to diverge in the vicinity of a charged current loop, and it is remarkable 
that the electrostatic potential plotted in Fig. [2] appears to be finite. This turns out to be a general feature of 
the zero-EMF solution, where both the electrostatic potential and the electric field remain finite at the loop in the 
Boyer-Lindquist coordinate system. As a consequence, loops of this type are useful when combined into an equatorial 
disc of constant potential. 

Consider the leading terms in the multipole expansion ( 18 ) which give a divergence as r approaches rp in the 
equatorial plane, 

^ - 2^^'"^7n E 'f7^Muo)Pi,r{uo)Pm'-' - E \^Mi{u,)P,{u,)P,\Ofx' (49) 

ri^(r+ - Af) ^ /(? + 1) ri^(r+ - M) ^ Z(/ + 1) 

where / = 47r/a — 2q, g ~ 47r/(ro + a?) — 2aq and a; is a dummy variable which has a; — >■ 1 as r — >■ ro- Standard 
formulae for the large I limits of the Legendre functions give. 



Pi 



MQiiuo) ^ -i-(r+-M)A(ro)-i/2 (50) 



2/ + 1 



Pi,rMQi,rM - ^|^(r+-Af)A(ro)-3/2 (51) 



The divergent terms in the sum are then 



A ^ V P,(0?.^ - V (52) 



4 ^ ^ 2A(ro)V2r2^Z(; + l) 
The sums can be expressed in terms of hypergeometric functions 



£PK0)V = F(l/2,l/2;l;x2) (53) 



1=0 
oo 



E l(7Ti) ^ ^' "^'^ ~ 2 ^' ^^^^ 

The hypergeometric function F(l/2, 1/2; 1; a;^) diverges as a; — 1. The potential is convergent if, and only if, the 
coefficient of this hypergeometric function vanishes, and that happens when 

AirlMa , , 

ro — 2m 



Consequently, in view of Eq. (45), the zero-EMF solution with vanishing total charge has a finite potential. 



8 



In an inertial frame, Gauss' law implies that the electric field of a charged loop diverges, and so a loop can only 
have a finite potential if there exists an inertial frame in which the charge of the loop vanishes. The appropriate frame 
is the one with orthogonal basis vectors 



e/ = (1,0,0,0) (56) 
,0,0,l) (57) 



—2Mar svc? i 
E - 2Mr 

This is the stationary frame, rather than the usual co-rotating frame of the Kerr black hole. We conclude that the 
zero-EMF solution with vanishing total charge also has vanishing local charge in the stationary frame. 



V. SPECIAL LIMITS 



The charged ring solution is likely to be most useful for giving the electromagnitic fields in vacuum regions close to 
the hole. In these regions a dipole truncation of the electromagnetic field can be a used as an approximation to the 
full field. The solution is especially useful for locating the innermost stable circular orbits [4: and for examining the 
production of energetic particles by the Penrose process. 

The effective ergosphere in which the Penrose process can take place extends far out along the equatorial plane 
with a relatively narrow height along the rotation axis [55]. Therefore, the fields in the equatorial plane can be used 
for studies of both stable orbits and the Penrose process. For a current loop around an extreme Kerr solution it is 
possible to obtain the vector potential in the equatorial plane in a useful closed form given below. 



A. Dipole and quadrupole truncations 

Starting with the inner solution, the Coulomb and Dipole terms combine into 



At = sin^ e+{Q- 2aMB) - ^) (58) 

BA clt 
A^ = —sm'^9 + iQ-2aMB)—sin^e (59) 



where A — {r^ + a'^)'^ — Aa^ sin^ 9 and 



The -B-dependent terms in each field component represent the fields of a black hole immersed in a constant magnetic 
field of strenght B Fig. [2] shows the dipole approximation for the zero-EMF solution with an extremal hole and 
a loop of radius 6M. The approximation agrees remarkably well for the At component and reasonably well for A^. 

Dipole truncations of the inner and outer general solutions for the vector field do not match smoothly on the sphere 
at r = tq. However, a combination of the inner and outer dipole truncations of the loop solution will at least give 
a reasonable approximation away from r = tq. For large r, with zero net charge the leading terms in the vector 
potential are a combination of a magnetic dipole and a charge quadrupole. 



At ~ ^P2{cos9) (61) 



A^ sin^6l (62) 



The dipole and quadrupole strenghts are 



-^{r+~Mfl3l d^'^{r+-Mfp;-^a{r+-Mff3l. (63) 



For the zero-EMF solution, these are given by 



nlroA^ (64) 
ro -2M ^ ' 

ro-2M 
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For large loops with rp >> r+, these reduce to the ordinary magnetic dipole strength m = Itttq and an electric 
quadropole moment d = —TrlaMro- 

B. Extreme Kerr limit in the equatorial plane 

In the extremal Kerr limit the parameter m — ^ oo and the asymptotic forms of the Legendre functions in the general 
loop solution can be used. The mode summations can be done in a similar way to those of Sect. |IVB[ and the results 
expressed in terms of hypergeometric functions, 

^1 1 / r-M ^ 



/ M\ ,,,,,, . s / ^^ 2ttIM^ 1 Q 

A= 1-— (C(ro)i^i(r) + ??(ro)i^2(r)) + (67) 



In the equatorial plane for r < tq 

A- ,„ 

r J rQ — 2M r r 

where 

r( \ 27r/Af 2 + g(ro ~ 2M) 

C(-o) = > (68) 

2nIM\ro + M)-qNP 
^ 2.o(ro-M)^ • ^''^ 

The outer solution for r > is given by the same expression provided that values of the hypergeometric function 
on the branch cut are defined by taking the average of values just above and below the branch cut. The zero-EMF 



solutions with vanishing total charge can be obtained by fixing the black hole charge to Q = —q and using Eq. (45 1 



^lAPF^jr) ( M\ 2i,IM^ 1 
* (ro - 2M)(ro - M) \ r J ro~2Mr' ^ ' 

This is finite at the ring, as we discussed in the previous section. The simple form of the solution accounts for the 
high accuracy of the dipole approximation. 

The axial part is not quite so simple. The axial vector potential is 

A<j, = ,7^^ 1 F2[r) + + 2ttIM (71) 

where 

Ci (r) = 2 TT Iro (ro - M) + (2 tt IM^ - 2 tt /ro^ + 4 ttIMto - qM) Mr^ 

+2 TT Iro^ (ro - M) Mr - inlroM* - q{ro - 2M)roM^ (72) 
C2(r) = 2n Iro {ro- I) r'^ + {2tt IM^ -2tt Ira^ + in IMro- qM) Mr + 277 1{ro + M)M^ - qM^ (73) 

C. Extreme Kerr with disks 

The availability of a closed form solution for equatorial current loops around extreme Kerr black holes allows us 
to examine what happens with thin equatorial disks. To construct models of this general type we must impose some 
physical condition on the disc. A suitable choice would be for the electric field to vanish is some specially chosen 
rotating frame with angular speed a;(r), 

At,r + i^A^^r = (74) 

This allows us to solve for the charge density q{r) given a current distribution I{r). as an example, we shall examine 
the simplest case, a; = 0, and find the charge distribution in a disk which has zero potential in the Boyer-Lindquist 
frame. 
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We can obtain the solution inside the disc by combining the solutions (67) for rings, 

A - f 1 - -) r [C(/)^^iM + v{r')F,{r)] dr' + ^ H ^"^dr' ~ ^, (75) 
\ r J J^^ r J^^ To - 2M r 

where ri is the inner edge of the ring. A zero potential solution can be obtained by having the first integral vanish 
and by fixing the black hole's charge 

Q-r^A'. ,70) 

In order to find out when the first integral vanishes, we use the identity AFi = (XF2)' + 2F2 to eliminate Fi. After 
integration by parts, this leads to two conditions, 

C(ri) = 0, 4C(r) +477(r) + (r-M)C(r)^^ = 0. (77) 

These can be regarded as equations which determine the charge density for a given current density, by substituting 
Eqs. (68) and (68). The first condition is the finite-potential condition. 

The second condition is a first order differential equation, 

[r{r - 2Mfq] ^ + M'^{r - 2M)(47rr/),^ = 0. (79) 

The solution is 



This is only valid for ri > 2M . Unfortunately, this appears to be the only case of ( 74 1 which seems solvable analytically. 



VI. CONCLUSION 



The main result of this paper has been to correct the published solutions to Maxwell's equations for the electro- 
magnetic vector potential near to a Kerr black hole surrounded by a charged current loop. The vector potential is 
a necesary ingredient for the determination of charged particle orbits close to a black hole immersed in a magnetic 
field. These orbits can be used to determine the inner edge of an accretion disc. They could also be used to extend 
recent work on energy flux generation using the Penrose pair creation 77 — >■ e"'"e^ and Penrose Compton scattering 
[29] to include magnetic fields. 

Current loops can be superposed to form non- vacuum solutions |16) . In these situations the loops will be charged. 
For example, an equatorial conducting disk with angular speed ujd would have a zero charge density in the disk frame, 
but a non-zero charge density in the Boyer-Lindquist coordinate frame. This can be modeled by a superposition of 
charged current loops in a similar way to the analysis in Sect |V C| It may also be of interest to combine the fields due 
to a current along the rotation axis with the current rings to produce a toy model which has a current circulating out 
along the equatorial plane and back down the rotation axis. 
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